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Abstract 



This paper is devoted to the hnear stabihty analysis of the whole spectrum of 
static hedgehog solutions of the Skyrme model on the three-sphere of radius 
L. These solutions are described by profiles F{;ip) that satisfy the equation 
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where -0 G [0,7r], -F(O) = 0, -F(vr) = Qvr and Q is an integer which has 
the interpretation of topological charge. This work is the continuation of 
the paper j3] where these solutions were classified. It turns out that only 
solutions that in the limit L —>■ oo tend to skyrmions (localized at the poles) 
are linearly stable. The other solutions are unstable and, for a given solution, 
the number of instabilities, for L sufficiently large, is equal to the index of 
a harmonic map to which this solution tends pointwise in the limit L — > oo. 
Such solutions, which in addition have a definite parity, undergo a transition 
by +1 in the number of instabilities as L grows. Due to the instability, 
new solutions, with spontaneously broken refiection symmetry, are born by 
bifurcations. In the case of the 1-skyrmion this critical phenomenon can 
be fully described analytically. This allows of prediction the unique series 
expansions for the profile of the 1-skyrmion and for its energy, though their 
expansion coefficients are not given in a general form, and read respectively 
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The above expansions are valid in some right neighbourhood of the critical 
radius L = a/2 at which the 1-skyrmion (with broken refiection symmetry) 
is born by a 'costing no energy' excitation of the marginal stability mode of 
parity +1 on the identity solution, which has parity —1. To the author's best 
knowledge the series were not given in literature so far. A similar mecha- 
nism of spontaneous breaking of parity is also observed when other solutions 
appear by bifurcations from symmetric solutions. 
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1 The Skyrme model on 

The Skyrme model is defined by the action which written in generally co- 
variant form reads 



where U is the basic chiral SU{2)-vaAued scalar field and K^^ = iU^dfJJ 
is a Lie algebra valued four-current. We use the convention (+,—,—,—) 
for the signature of the canonical quadratic form associated with a metric 
tensor. The first integrand is the cr-term, the second one was introduced by 
Skyrme in [H] to ensure the existence of solitons, but its shape may be also 
determined by using geometrical arguments [Zl,j3]. Both these terms have 
mutually inverse scale dependence, thus compete against each other. Unlike 
in the Skyrme model on flat space, in addition to the natural size of a soliton 
(ef^)"^, we have another natural size at our disposal which is introduced by 
the radius R of the base three-sphere. Their quotient L = ef^i? (which is 
a dimensionless number as the vr) must have nontrivial consequences to our 
model since the parameter L, as a remnant of the mentioned competition, 
enters the equations of motion. 

In j3] we found and classified all static, finite energy, spherically symmet- 
ric and equivariant mappings F from the base space - the three-sphere of 
radius R, to the target space - the group SU{2). These mappings are regular 
solutions of the equation 



which is singular on boundaries. We utilized the canonical correspondence 
between points (\E', 0, $) on and elements of SU{2) (the last is metrically 
equivalent to S^), that is if t/ G SU{2) then 
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[U] = I V^d'x(^^f'^TT{K,Kn + ^TT{[K,,K][K^,K'^])y (1.1) 
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The base three sphere is endowed with the standard spherical coordinates 
'(9, (fi) in which the line element reads 



ds^ = (dV-^ + sin^ (dt?^ + sin^ dd^'^)) , e [0, tt] x [0, tt] x [0, 27r] 

The equivariance of the mapping means that 

F: 53 9 ^ = F(V'),e = t?,<l' = (^) G 5^ (1.3) 

thus is related to the Skyrme's ansatz, which is popularly called the 'hedge- 
hog' ansatz. 

The number of solutions of equation ()1.2|) increases as L grows. In the 
limit L — > 00, a given solution tends to a configuration which is composed of 
two multi-skyrmions S„ and S„, respectively, localized at the north and at 
the south pole of the base three-sphere, and a 'harmonic map' in between. 
These solutions are denoted symbohcally by S„HpS„. This means that each 
solution tends pointwise to a harmonic map hp of Bizoh IJ in some interval 
6 < ip < TT — e, and e —>■ a.s L —>■ 00. If, in the vicinity of the north pole, 
one introduces the radial variable r defined by 1 ^ ^/^ oc rL~^, then F[r) 
approximately satisfies the equation of the Skyrme model on flat space [S]. 
In the vicinity of the south pole one can proceed in an analogous way. 

In what follows, for convenience, we will often refer to F'{0) and -F'(vr) 
as shooting parameters and also to be in connection with 3j where this 
nomenclature was introduced due to the specific method used for finding 
solutions. 

The solutions were examined to see which static configurations of the 
field F are possible at a given size of the base three-sphere. Energies of these 
configurations are given by the functional 



U[F]= w{ij)U{ij,F{ij),F'{ij))di:, w (tP) = in sm'^ (1.4) 
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As the unit of energy we chose ef^^/2. Due to the spherical symmetry 
imposed on solutions, the volume element on was integrated over S"^ giving 
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rise to the weight function w{ip) = Air sin^ ip . This function is important for 
further considerations since it is natural to use it as the weight function in 
the vector space of spherically symmetric functions defined on S^. 

2 Linear stability analysis. Some general in- 
troductory remarks 

Spherically symmetric, equivariant equilibria of the functional (jl.lj) . that is 
(time independent) critical points of energy functional ()1.4|) . may be stable 
or not in dependence on whether they minimize the functional l^[F], which 
can be interpreted as potential energy. To check stability of a map F which 
is a solution of equation 6U[F] = we compare its energy U[F] with en- 
ergy of its deformation U[F + e^], where ^ is a trial function satisfying some 
conditions (to be stipulated later) and e is a small real number, which is 
used as the parameter of the formal Taylor series expansion of the functional 
U[F + eC,]- In this paper we are restricting ourselves to examining only the 
behaviour of the term 0{e'^) of the expansion, by which it is tantamount to 
the linear stability analysis of solutions F (this procedure is analogous to ex- 
amining local extrema of ordinary functions for which /" 7^ at extremum). 
Our aim is to explain bifurcations of solutions in our model and it suffices 
to examine linear stability alone. This is also the reason we ignore in this 
paper the cumbersome problem of marginal stability when a perturbation 
does not change the term O(e^). In what follows we will often skip the ad- 
jective 'linear' for brevity. For an equilibrium solution to be stable we thus 
require that in its vicinity the second variation of V([F] be positive definite. 
The positive definiteness of 6'^U[F] is necessary and sufficient condition that 
a critical solution was stable in the domain of spherically symmetric pertur- 
bations. Nevertheless, it should be kept in mind that a solution may turn 
out to be unstable if the constraining conditions imposed on symmetries of 
perturbations are weakened. 

Below a system will be said to be energetically stable if it satisfies the 
above mentioned 'abstract' criterion of stability. This is to distinguish it 
from the intuitively well understood notion of dynamical stability. The last 
is being associated with the evolution of a system in time. This evolution is 
determined by a kinetic term of the total energy functional of the system. 
Dynamical stability of a system in equilibrium is understood as remaining 
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arbitrary close to the equilibrium state if only the perturbation is sufficiently 
small. The measure of the perturbation may be determined by the amount 
of energy the perturbation introduces into the system. 

In appendix El both these notions of stability are compared one with an- 
other. This comparison was undertaken to see the differences and whether 
they can affect the general statement that a solution is linearly stable or 
not. To some extent these notions turn out to be equivalent but some theo- 
rems on spectral analysis are required when more than one mode of instability 
appears, since the two approaches result in different Hilbert spaces of pertur- 
bations (the difference comes from weight functions). This will be elucidated 
in the following paragraphs. 

2.1 Formulation of the problem 

From now on one can forget about the genesis of energy functional ()1.4|) and 
consider it as defining an elliptic problem for a spherically symmetric scalar 
field F 'living' on the unit three-sphere S^. Then L plays the role of a dimen- 
sionless free parameter. Classical solutions of the model, that is the critical 
points of functional ()1.4|) . satisfy equation ()1.2|) . This equation is necessary 
but not sufficient condition for a solution F to be a local minimum of the 
functional. To see whether it is a minimum it suffices to check if for a class 
of perturbations the energy U[F + e^] is greater than U[F]. The perturba- 
tion ^ can be normalized and e denotes its (small in the sense of some scalar 
product) amplitude. If the energy is increased for all perturbations within 
some given class we say that solution F is energetically stable with respect 
to this class. For our purpose, among other things, we have to assume this 
class to be composed of all perturbations which do not move the considered 
fields F away from the topological sector and do not affect their spherical 
symmetry. The condition of symmetry for perturbations allows us to reduce 
our problem and to treat it as a one dimensional since one can consider the 
differential equation ()1.2|) as such, regardless of its origin. Thus, one can 
reformulate this problem as follows: in the domain G [0, vr] find all regular 
solutions of (|1.2p and check their stability. To do this one may look for some 
variational principle which reproduces the equation. But this may be done 
in many ways, e.g. by inclusion of more spatial dimensions and time. 

To give an example we can treat the integral ()1.4|) as a generalization of 
the energy functional of harmonic mappings between three-spheres. Then 
to check energetical stability of these generalized mappings we would follow 



4 



PP by choosing simply 4tt sin^ ip as the natural choice of the weight function 
in space of spherically symmetric perturbations on S^. On the other hand, 
since originally we obtained p.2|) from the Skyrme model on S^, we would 
write down the wave equation for spherically symmetric, time dependent 
and equivariant solutions F{ip, t). Static solutions of this wave equation, the 
equilibria, would also be solutions of (jl.2p . We would then follow the stabil- 
ity analysis made for spherically symmetric static solutions of the Einstein- 
Skyrme model (of self-gravitating skyrmions) [2\. This would correspond to 
the stability analysis in dynamical sense. Both the approaches would lead 
to the same Sturm-Liouville operator acting in the same space of admissible 
functions of perturbations but endowed with different scalar products - thus 
we would obtain different Hilbert spaces. The difference would come from 
the weight functions - the first method would give simply w{ip) while the sec- 
ond would give a complicated solution-dependent weight function. To make 
things worse, even in the case of one space dimension one can relate solutions 
of ()1.2|) to equilibria of many 1 + 1 dimensional field theories that differ one 
from another in kinetic terms. In fact, to examine energetical stability of a 
field -F(V') in one space dimension one has to introduce some weight function 
for which, among other things, the most natural choice is the standard vol- 
ume element on the interval [0, vr]. To this weight function there corresponds 
some 1+1 field theory for which F{iIj) is an equilibrium solution, and whose 
kinetic term is such, that the equation for eigenvibrations around equilibria 
is the same as the Sturm-Liouville equation which originates from the en- 
ergetical stability analysis. In this case both the dynamical and energetical 
stability analysis lead to the same results (to see the correspondence between 
the weight function and the kinetic term the interested reader is referred to 
appendix ini). 

Consequently, it arises the problem of a definition of an appropriate uni- 
tary space of admissible perturbations, that is the class of functions which 
should be used to vary the functional U, and of the choice of a weight func- 
tion used to define scalar product in this space. From the above analysis 
it follows the arbitrariness in the choice of weight functions. The conse- 
quence is the question how this arbitrariness affects the spectrum and qual- 
itative results of stability analysis. Can we simplify the analysis of ener- 
getical stability of hedgehogs by choosing simply g = wl (This choice of 
the weight function is due to the 'simplicity and naturalness'). Put dif- 
ferently, we would like to know whether the qualitative predictions of this 
analysis are norm-invariant and, in particular, if the dynamical stability anal- 
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ysis of time-dependent fiedgeliogs (wliicfi produces more complicated Sturm- 
Liouville equation) would give the same qualitative results? 

As an aside, we remark also that from the point of view of the functional 
(jl.ll) the condition of spherical symmetry for perturbations may turn out 
to be too stringent. In fact, a solution which is stable under the stability 
criterion we assumed, might turn out to be unstable in wider sense since, 
due to some infinitesimally small nonspherical perturbation, it would evolve 
to another configuration. With this reservation in mind we assume the con- 
straints of spherical symmetry on admissible perturbations and, in addition 
to vanishing on the poles of the base S^, we assume them all to be pointwise 
continuously differentiable. The last condition follows from the requirement 
of continuity of the energy functional (jl.4|) . 



2.2 The Hessian and its weight function dependent 
spectrum 

The Hessian measures energies of excitations ^ around a classical solution F 
and is quadratic with respect to perturbations ^. To avoid confusion between 
different conventions, we define the Hessian as the coefficient 6'^U[F]{^,^) in 
the series U[F] + SU[F]{^)e + S'^U[F]{^, ^)e'^ + ... which is the formal series 
expansion of the functional U[F + e^] with respect to the variable e, hence 



(2.1) 



where 



L + 



2 sin^ F 
L sin^ i/j 



4 sin 2F 
L sin^ ip 



+ 



sin^ F 2 cos 2F 
2cos2F)-^ + — ^ 
sm ip sm ip 
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This is tantamount to defining an operator Ti which acts in a linear space of 
admissible functions A whose elements are used as perturbations. If |0 ^ 
then 7i is defined according to the formula {^\7i\^) := 5^W[F](^,^). For the 
form 5'^lA[F]{S,,i) to be continuous we require A to be composed of piecewise 
continuously differentiable functions. Suppose that we have succeeded in 
finding some countable and complete set of eigenvectors G ^ to the 
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corresponding eigenvalues Aj of the operator 7i. To ensure this we require 
the form (?7|?-^|^) to be symmetric and the function space A to be endowed 
with some scalar product which, for any two vectors |.^) and \ri) from 

A, is defined by the integral 

TT 



The weight function g has to be positive for ip G (0, vr) and fulfill some other 
conditions to ensure the existence of the above integral for all elements of A. 
Now each perturbation C, & A can be uniquely decomposed in A and thus 
written equivalently in the form af an infinite series \0 = - Ci\^i). To excite 
the perturbation the amount of energy {^\T-C\0 = Xli ■^j'^i ll'^^llg required. 
The representation of the function ^ in the space A is unique and given by 
Ci = ('CilS'lO ll^illg ^- Thus, to decide if a solution F is stable it suffices to 
check the positive definiteness of the quadratic form (^IT^IO; to check 
if all eigenvalues are positive. From this it results the following criterion of 
linear (energetical) stability: 

Criterion of stability 2.2.1 For a solution of il.^) to be linearly stable 
in the domain of spherically symmetric, pointwise continuously differentiable 
and vanishing on boundaries perturbations, it is necessary and sufficient that 
the resulting spectrum of eigenvalues of Hessian \2.1]) . that is of the second 
variation of functional \1.4^ , evaluated at this solution, was positive. 

Now the problem of a choice of the appropriate scalar product arises. Any- 
way, before deciding this, we will find the operator H and carry out its 
spectral decomposition assuming, temporarily, an arbitrary scalar product. 
We assume that all eigenvalues are enumerated in such a way they form a 
nondecreasing sequence {Aj}. Then it follows, analogously as in the case of 
finite dimensional quadratic forms, that each element A, may be character- 
ized by a process of consecutive minimizations |3]. The lowest eigenvalue 
Ao is defined as the global minimum of the functional • \\0\g'^ 

the domain A. By we denote the minimizing ^. The n'th eigenvalue is 
a minimum of the functional ■ ||^||~^ under the assumption that 

iOdl^k) = 0, k < n and the minimum is attained by ^ = Applying this 
to the functional 

Af£l-=i^ eeA 
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it is straightforward that, at least for the global minimum ^ = ^o; the equation 



O.MKl.^WMM), w„e.e A:^AK1 (2.2, 

ml 

and consequently the equation 

or equivalently the equation 

1 6{wT) 



must necessarily hold since > 0. The foregoing functional derivative 
involves F'^ip) which, for solutions, can be expressed by F and F' using 
fra . hence 

-^(^Tii Lsin^^ + |sin2F(V^)^ (2.4) 



dip \ \ L J dip J 

L sm ip + j; sm 

where 

A(V^) = (2cos2F(^) - l)sm'^ F{ip) 
\ sm ip ) 

(F'(ih))'^ 

+F'{ilj)sm2ilJsm2F{ip) + ^ ^^^^ (l - cos2F(^) (l + L^sin^V')) . 

Equation ()2.4|) together with the class of admissible functions A satisfy 
the general requirements for the Sturm-Liouville problem that the theorems 
proven in jl] would be successfully utilized. It follows that the condition 
i^lolCk) = 0, k < n in finding the other consecutive minima A„ is sufficient 
that equation ()2.2p would hold in general. It also follows from ^ that the 
respective minimizing C,n exist and are the same as solutions of equation ()2.4|] 
which is the necessary condition that the first variation of ()2.1|) vanished. It 
also follows that the corresponding eigenvalues A of (j2.4p are the same as the 
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minima A^. Moreover {^i\T-C\^j) = KWdWgSij and the denumerably infinite 
and, by construction, mutually orthogonal set of eigenfunctions is com- 
plete in A. Using the metric form we normalize these functions to the 
unity. Thus, we have managed to construct the Hilbert space {A, \\-\\g) in 
which form an orthonormal and complete base and in which the operator 
H is consequently represented by the formula 

oo 
A:=0 

The second order linear self-adjoint expression on the left in (j2.4j) is just 
the element {ijj\T-C\C,) ■ As an aside, we remark that in this way we were led 
straightforwardly to the differential operator 7i which is naturally born as 
self-adjoint. Due to the boundary conditions imposed on elements of A this 
operator is also hermitian. 

Before we utilize the foregoing formalism to carry out the stability analy- 
sis of solutions of equation p.2|) . we examine the problem of deciding which 
weight function should be substituted in place of g. It is straightforward for 
spherically symmetric functions defined on the unit three-sphere to normal- 
ize them using the natural volume element w{ip), i.e. to substitute g = w 
into (j2.4j) . In this paper we follow this natural choice, especially, as then the 
comparison with the results of harmonic maps between the three-spheres 
may be done. On the other hand, one would equally well argue for another 
choice. In fact, the functional ()1.4|) defines a problem of finding extrema in 
the function space composed of functions of one variable ip ^ [0,7r], there- 
fore the natural choice of the weight function would be simply g{ip) = 1. 
Another possibility is given by the following physical argumentation. The 
weight function should be chosen in such a way that the resulting spectrum 
of the Hessian could be directly interpreted in terms of frequencies of eigen- 
vibrations of our system and, as such, should be determined by the kinetic 
term alone. The appropriate g was constructed in appendix lEl 

It seems there is no sufficiently strong criterion which would determine 
the appropriate weight function. This signalizes that it is rather quite ar- 
bitrary which weight function should be chosen to normalize the function 
space of admissible perturbations, as long as the axioms of scalar product 
are satisfied. The situation with the arbitrariness of the choice of 'g' re- 
sembles a sort of gauge freedom. However a nontrivial change in the scalar 
product inevitably affects the spectrum and the respective eigenfunctions. 
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Consequently, this changes the Hilbert space, so that 'observables' can not 
be unitary transformed one to another. Nonetheless, it is very plausible this 
is not a real obstacle and may be successfully cured. Some arguments it is 
really the case are given below. 

From the minimizing properties of eigenvalues it trivially follows the ob- 
servation that the positive definiteness of the Hessian is universal, i.e does 
not depend on the specific weight function assumed. It is also clear that the 
sign of the lowest eigenvalue is universal as well. Thus, to answer the general 
question if a solution is stable, it is quite arbitrary in which scalar product 
the function space A is endowed. This is in agreement with the intuition that 
stability is not of dynamical origin. The problem arises when one wants to 
know how many instabilities a particular solution possesses. The question is 
whether the number of instabilities, i.e. if the number of negative eigenvalues 
(which is always finite) is norm invariant? Actually, a simple argument that 
the number of instability modes is norm-invariant can be constructed. 

It turns out that at a critical L^, at which a new solution of equation 
(ll.2j) appears or bifurcates from an already existing one, the solution has 
an eigenvalue which vanishes. Let denote it by A„ where n is the number 
of nodes of the corresponding eigenfunction Since A„ = it is clear 
that this eigenfunction is universal at L = Lc, i.e. does not depend on a 
specific weight function g (to preserve the asymptotics of solutions of ()2.4p 
we may assume that weight functions vanish on boundaries). Moreover, the 
number of nodes of which is also the number of negative eigenvalues, 
can not be affected by a continuous change of any coefficient of the Sturm- 
Liouville equation ()2.4|) . Due to continuous dependence of eigenvalues on the 
coefficients a negative eigenvalue must remain negative for all L such that 
Lc < L < Lc+e where e is positive and sufficiently small. From the continuity 
and universality of Xn{Lc) it also follows that the sign of A„(L) must be 
universal for < L — Lc < e. Otherwise there would exist such two weight 
functions ^'(0) and g{l) for which, respectively, A„(0, L) < and A„(l, L) > 0. 
Then for some a G (0, 1) the equation ()2.4|) with g = ag{l) + (1 — (y)g{0) 
would have vanishing eigenvalue A„(a,L) = which, in turn, would not be 
universal, a contradiction. 
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3 Linear stability analysis of solutions S^H^S^. 
Critical phenomena and spontaneously bro- 
ken reflection symmetry 

To summarize, we have reduced the problem of examining the hnear stabihty 
of solutions of equation ()1.2|) to the Sturm-Liouville problem for the differen- 
tial operator ()2.4|) with the condition that its eigenf unctions ^ [ip) must vanish 
on boundaries ^(0) = = ^(vr). As thea weight function we chose simply 
the volume element on - g{ip) = 47rsin^'?/' - as we gave some arguments 
that qualitative results of our analysis can not be affected by this special 
choice. In this way we also avoid the unnecessary complications introduced 
by the weight function implied by the kinetic term of the Skyrme model 
on S^. Due to global regularity of solutions of equation ()1.2p the boundary 
points are regular singular points of equation (|2.4p . Thus solutions of (|2.4|) 
in the vicinity of boundaries can be written as generalized series. From the 
secular equation of (|2.4p it follows that the first independent solution is a 
Taylor series for which ^ vanishes at a boundary point and thus can be made 
vanishing on the second boundary point to satisfy the requirements of ad- 
missibility. We also note the very useful for numerical integrations fact, that 
(if not degenerate) the eigenfunction corresponding to the eigenvalue 
has exactly k nodes within the interval (0,7r) (i^o has no nodes inside) and 
the sequence {Xk} is nondecreasing and divergent. 

3.1 Hk in the limit L ^ oo, harmonic maps 

Let F be the solution (then F' is finite for all ip G [0,7r]) and take the 
limit L — cxD in ()2.4j) . In this way we reproduce the differential equation for 
linear perturbations of harmonic maps h^ of Bizoh which was analyzed in P 
(the author used the conformal variable x = In (tan {iIj/2))) 



The equation can be solved analytically for the identity solution, the har- 
monic map hi (the case of the vacuum Hq is trivial) then 



iC' sin^ ipy + 2 cos {2F)C = M sin^ ^, A = lim 




1 



sin ip Cf {cos ip), 



uj] = ^2 + 4/ - 1 
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where (and respectively uj) denote the Vth (where / = 0, 1, 2, . . . ) eigen- 
function (eigenvalue) around the regular solution h^ and Cf are the Gegen- 
bauer polynomials. The map hj. has exactly k unstable (S'0(3)-symmetric) 
modes (for detailed information see P). 

3.2 Stability of the vacuum solution Hq 

It is clear that the solution Hq (Hq = kir, A; an integer) whose energy is zero 
must be stable, since energy integral ()1.4|) is bounded from below by zero. 
Hence any perturbation can only increase energy, therefore the vacua are 
stable. Nonetheless, we carry out the calculations since the spectrum will be 
used to interpret spectra of perturbations around S„ or S„S„ . This is because 
skyrmionic solutions tend pointwise (but not uniformly) to Hq in the limit 

L — > CXD. 

The substitutions F{ip) = kn and Vl — x'^u{x), where ip = 

arccos(— x), reduce equation (j2.4p to the eigenvalue problem 



which is the Gegenbauer equation for u. Due to the boundary conditions 
imposed on ^ the function u must be bounded everywhere in the closed 
interval [-1,1]. This is possible if A = A; = L{f + 4/ + 3), / = 0, 1, . . . . Thus 
the eigenvalue problem for Hq is solved and given by 



where Ci are the Gegenbauer polynomials. Written explicitly, the few first 
normalized modes of the vacuum reads respectively 






sin2?/' (1 + 4cos2?/') . 



12 



3.3 On the instability of which gives rise to the 
appearance of Si by a critical phenomenon at the 
critical radius L = \/2 

Substituting into ()2.4|) F{ip) = ip we obtain the eigenvalue problem 

-(sin2^e'(^))' + 2 (^1 - 2^±^sm'ij^ m = -A-sin^^A^W (3.1) 
which is solved analogously as before resulting with the Gegenbauer equation 
(1 - x^)u"{x) - 5xu{x) + { j^2^2 ^ + " ° 

whose solutions are bounded everywhere in the closed interval [-1,1] only if 
A = A; where 

Xi = yif + Al + 1) + L{1^ + 41 - 1). 

The respective eigenf unctions of ()3.1|1 are the same as for Hq. (If we had used 
the weight function resulting from the kinetic term of the Skyrme model we 
would have got = LXi/{L'^ + 2) with A; as above, see appendix IB]). In 
the exceptional cases of Hq and the respective eigenvectors do not change 
with L for the weight function assumed. It is no longer true in generic case 
of weight function. (The other solutions of (jl.2p . known only numerically, 
change continuously as L increases, thus so do the eigenfunctions). 

All A; with I > 1 are positive for every L. The map Hi is stable for 
L < \/2 and consequently a minimum of the functional ()1.4|) . For L = 1 it 
also saturates the Bogomolnyi bound E = 127r^. For L > \pl the map Hi is 
no longer stable and by this the critical value L = is distinguished. When 
the threshold of stability is passed, the new solution Si appears. Due to the 
reflection symmetry of equation ()1.2j) : F{iIj) it — F{tt — ip) it also appears 
the 1-skyrmion Sf which is localized at the south pole for large L (in what 
follows we will be using this notation for that 1-skyrmion). Note that Si is 
not reflection symmetric unlike Hi. Both Si and Sf bifurcate smoothly from 
Hi and exist for L > \/2 tending continuously to skyrmions localized at the 
poles of the base three-sphere in the limit L oo. In figure ^ the evolution 
of spectra of eigenvalues of ()2.4|) for Si and Hi are shown. Due to positive 
definiteness of the Hessian at F = Si it follows that the skyrmionic solution 
is stable (it possesses also the lowest energy). 



13 



Figure 1: Evolution of eigenvalues for Hj ( -) and Si ( ). We use 

the definition := 2L~^ where L is the radius of the base three-sphere, and 
instead of eigenvalues A we show the rescaled values \/L that are finite in 
the limit — > 0. At the critical = 1 {L = a/2) the lowest eigenvalue of Hi 
becomes positive, thus Hi is stable for > 1 and unstable for < 1. Due 
to marginal stability of Hi at the critical = 1 the skyrmionic solutions Si 
and Sf appear bifurcating from Hi. The skyrmionic solutions are stable for 
<1 and do not exist for > 1. 



In what follows, we explain the numerical observations. In a sufficiently 
small right neighbourhood of the critical L it is energetically preferable to 
excite the mode on Hi since then the energy is diminished. It is also the 
only mode with negative energy at our disposal. The exceptional role of the 
mode ^0 is also reflected in the fact that it is the invariant solution of the 
eigenvalue problem ()2.4|1 at L = \/2 with respect to a change of the weight 
function. Thus the solution has to play a nontrivial role in the vicinity of the 
critical L. Close to the critical value L = ^J2 the plot of shooting parameters 
of the map Si (figure |21) is (at confidence of 0.95) perfectly described by the 
fitting curve 

L - = 5 + c(a - ao)^(l - Ci(a - a^) + c^ia - 1)^ + . . . ), (3.2) 
Oo = 0.9999997 ± 4 x 10~^ 6 = 2, 0000007 ± 0.0005, 
c = 0.2592695 ± 5 x 10"^ 5 = -1.54 x 10"^ ± 10"^ 
ci = 0.6001 ±0.02, C2 = 0.50175 ±0.1. 
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Figure 2: Shooting parameters - a - (on the left) and energies - E = U[F] 

- (on the right) for ( ) and Si ( ) as functions of and L, 

respectively. We use the definition := 2L^^ where L is the radius of the 
base three-sphere. At = 1 (L = a/2) the skyrmionic solution Si bifurcates 
from Hi and exist for < 1. The curves of shooting parameters of Si and 
of Hi (for Hi a{L) = 1) cross vertically each other at the bifurcation point 
= 1 while the curves of energies (on the right) smoothly approach one to 
another. This characteristic singular 'fork-like' structure in both diagrams 
signalizes a critical phenomenon. 



This is the numerical confirmation of the hypothesis that in the vicinity of 
the critical L the shooting parameter a of the map Si possesses the fractional 
exponent behaviour 

|a- 1| oc (L- V2)^/^ (3.3) 

characteristic for critical phenomena. This behaviour is remarkable and it 
is desirable that its analytical explanation was given. In what follows we 
will show that this result may be reproduced analytically together with the 
calculation of the proportionality factor in ()3.3|) . 

In the vicinity of critical L = \/2 we can find three distinct formal Fourier 
series expansions of the function F{ip) — ip which depend on some small real 
parameter a 

oo 

F{il)) — ip a'^amia) sin {mip), 

m=l 
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where amO: are functions to be determined and a ^ as L ^ v^"*"- Its 
shape guarantees that as a — > then F^ip) — > Hi. The other solutions S„S„ 
or S„HiS„, which exist at L < a/2, are unfortunately too remote from the 
analytically known Hi to be examined in this way. The coefficients am(«) are 
defined by the requirement that F{ip) should satisfy equation ()1.2|1 expanded 
in a. These coefficients may be found by a recurrence process if we assume 
that they may be approximated by partial sums of the form 

N 

a(^)(a) = a„,fea2^ with := L - V2 > 0. 

k=0 

Since the eigenvalues of Hessian at Hi read ^2n(V') oc Ylk=o('^^~^^) ((^^ + ^)'^) 
and ^2n+i('^) X]fc=i 2^ sin {2kilj) it is clear that there exist formal (infinite) 
polynomials bm{a) of the variable with 6m (0) 7^ such that the above 
Fourier series may be recast into the form 

00 

m=l 

Before we proceed further, it should be stressed the fact that this procedure is 
complitely formal, since it presupposes a line of analytical properties that do 
not have to be valid. Therefore the procedure may lead to false conclusions, 
unless one proves its correctness. The first of the three series is simply zero 
and corresponds to the solution Hi. To the order of a^, to write only a few, 
the second is given by 




ou V ^ 29369 / r- 
F{i;)^^P+ I ^/_^-a2__^i65v^| asinV' 



1 ^ '^/^ 

sin2?/' + a^ sin3V' + o(a^) (3.4) 

22 88 

and the third formal series may be derived from the last by the refiection 

F{iIj) ^ tt — F{tt — ip). The reader is referred to appendix ^ to see the 

expansion up to a higher order. For L < \pl with o? = \pl — L, a similar 

procedure does not give any real series apart from the one which is identically 

zero and corresponds to Hi. Expressed in the base of the eigenvectors of the 

Hessian fl2.4p evaluated at Hi the foregoing series reads 



5\/2 9\/2 
F(V') ~ + 3vr^/ ■ <o(^) - ^vra2^i(^) + o{a'). 
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The question arises whether the formal series (j3.4j) has a nonzero radius of 
convergence in the variable a. It requires a rigourous proof. Here we give only 
a naive argument which makes it plausible that ()3.4j) is really the solution of 
()1.2j) (the reader is also referred to appendix^ to see how these coefficients 
behave). It was proved in jS] that finite energy solutions of equation ()1.2|) . 
in the vicinity of the singular point t/) = 0, are analytic functions of tl). Thus, 
every such solution can be written as an infinite power series whose general 
form reads 

The series contains only odd powers of (the series ()3.4|) as well) and the 
coefficients are rational (thus also analytic) functions of a and L. Substituting 
into the series L = a/2 + and a = F'(0) = 1 + ca + . . . calculated from 
(|3.4j) and next expanding in a we compare the result with (j3.4j) expanded in 
(the comparison was carried out up to tp-^^ and with positive result). 
Now we can compare the 'solution' ()3.4j) with the numerical results for Si. 

As one would expect this approach explains our numerical results. The 
formal series predicts in the vicinity of L > ^/2 the singular behaviour of 
shooting parameters of the map Si 



a 



- 1| = ■ (L - + o{{L - v^)V2. 



To compare with the numerically derived curve ()3.2|) we rewrite the above 
formula into the form L - = 0.2592725 ■ ■ ■ ■ {a - if + o{{a - if)). The 
comparison is also a sort of accuracy test of the numerical integrator used to 
solve equation (jl.2j) . In spite of the fact that close to the bifurcation points 
the numerics can not be very precise, this test gives quite good results. 

From the foregoing analysis it also follows that, to first order in the critical 
parameter a, the instability which causes Si (and Sf ) to appear by bifurcation 
from Hi at L = -\/2, is due to the mode whose eigenvalue becomes negative 
at L = -\/2- Thus the mode is in fact distinguished and this is due to the 
nonanalyticity of the function a{L) at the critical L = ^J2. The function sin-?/; 
was used as a trial function to verify the instability of the identity map in jH] 
where the Skyrme model on was proposed. Simultaneously, the authors 
expressed the belief that with a better trial function it should be possible 
to establish the instability of the identity map at L < It is clear that 
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this must be performed by using trial functions which are not spherically 
symmetric since, at is was shown before, is a local minimum of (jl.4p for 
L < \pi in the domain of spherically symmetric functions. 

The formal solution Si in the vicinity of critical L ()3.4|) can be utilized 
to find a similar series expansion of energy of Si. To carry this out up to 
we need to take into account contributions from the three first modes what 
gives (to see the expansion up to a higher order see appendix [XJ 

,,rn / M 2/^3^2 1 2 19^2 4 15627 

l/|S.(a)] = 12.' + -c? - —a' + —a" + 0(a') j . 

To compare with, we also expand W[Hi] 

W[Hi] = 12n^ + + - + 0{a')^ . 

Subtracting we get {U[S,{a)]-U[R,{a)]) / {12n'^) ^ -0.482118a^+0.491900a6+ 
0{a^) which is in agreement with the numerical result — 0.482(l)a'^+0.49(12)Q;^ 
Thus, the energy of Si behaves smoothly and its plot coalesce with this of 
Hi in a characteristic cusp at the critical point L = a/2 where both curves 
are tangent one to another (figure Ej). 

The approach presented here has the disadvantage that one can not be 
sure that the series ()3.4p is really a solution of ()1.2|) unless the appropriate 
proof is known (nevertheless the numerics was quantitatively and qualita- 
tively explained). Therefore we give qualitative and mathematically correct 
explanation (nonetheless quantitatively inaccurate). The evolution of the 
profile Si as L grows resembles very much the conformal deformation of the 
map Hi and it is well seen if instead of ip one uses the conformal variable 
X = lntan(?/;/2) [H]. It should be noted that the first time the conformal 
deformation of the identity solution was examined was by Manton in and 
it was in connection with the stability analysis of the identity solution, (our 
approach and motivation here is quite different). Here we compare the in- 
stability result with the conformal deformation of the identity to construct 
a function which would have a similar singular behaviour at the critical L. 
The conformally deformed Hi (denoted by Hf ) can also be decomposed in 
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the base of eigenf unctions of the Hessian at F = 



2 arctan (^e^ tan ^) = + ^ " + J P^o 
The energy ofHf can be found exactly and is given by 



^ V ^ cosh2(/5/2); 



To chose the conformal parameter as good as possible we can define it by 
the requirement that this energy, i.e. the function = W[Hf], was at 

minimum. For L < f'{(3) > thus f{/3) attains its minimal value at 
l3 = which is the energy of the map Hj. For L > a/2 the equation f'{(3) = 
is solved also for some P{L) ^ which corresponds to the energy 

Z^[Hf (^)] = Un' (V2-^^ < W[H,], L>V2. 

The function (3{L) has the analogous interpretation as a{L) — 1 for Si and 
also possesses the characteristic critical behaviour (3{L) ^ \/IL — V2) in 
the vicinity of critical L = ^/2. Thus again we have constructed a function 
which is singular at L = ^/2 and the corresponding profile F{tjj) (which is 
not solution of p.2|) ) has perfectly smooth energy. In the limit L —>■ oo 
this energy is finite and close to the energy of Si for which, in this limit, we 
have 127r^ ■ 1.23145. Moreover, in this limit -F'(O) of H.^^^^ has the analogous 
behaviour as the shooting parameter oq of Si, that is F'{0) ~ L. The shooting 
parameter of Si at the south pole behaves like for large L while for 
nf*^^^ it behaves like L~^. Anyway, in this limit, the solution Si is localized 
at the pole ip = thus Hf*-'^'', in a sense, well shows qualitative properties of 
Si and possesses the analogous behaviour at the critical point. This enables 
us to qualitatively understand the critical phenomenon. 

As an aside, we remark the fact, that the critical behaviour is observed 
neither in the Skyrme model on flat space nor in the model of harmonic 
maps between three-spheres when free parameters of these models - the 
characteristic size of a soliton in the first case, or the radius of the base 
three-sphere in the second case - are changed. This shows that by coupling 
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together two field theories we may produce new phenomena that are absent in 
the decoupled case. Dimensional parameters introduced by different models 
determine then some numbers whose value is usually crucial to the existence 
of different types of solutions. This fact is known when nonlinear fields are 
coupled to gravity |2]. 

In what follows, we shall say that F{iIj) has parity +1 or —1 if, respec- 
tively, 

/(0) = /(-0) or /(0) = -/(-0), 

where 

/(0):=F(vr/2 + 0)-^^^2^1±^, g [-7r/2, +7r/2]. 

The solution Si has indefinite parity, unlike Hi, from which it bifurcates. This 
phenomenon can be easily explained by analyzing, in the vicinity of L = -\/2, 
the decomposition of Si in the base of eigenfunctions of the Hessian evaluated 
at Hi. This is another remarkable feature of the critical phenomenon in the 
Skyrme model on S^. We postpone this explanation until the last paragraph 
where quite analogous behaviour is observed when other solutions appear by 
bifurcations. 

4 Linear stability analysis of solutions found 
numerically 

To classify solutions of equation ()1.2|) the notation S„HpS,„ was introduced 
fS] to stress the fact that in the limit L ^ cxo a given solution tends to a 
limiting configuration composed of a harmonic map h^, of Bizoh |lj localized in 
between the poles of the base three-sphere, to which the fiat space skyrmions 
are attached at the poles, respectively, the n-skyrmion is localized at the 
north pole {ip = 0) and the m-skyrmion is localized at the south pole. This 
configuration is characterized by a total topological charge Q which is equal 
to n + m or to n + m ± 1, respectively, if the solution contains an even or 
an odd harmonic map. The symbol S„HpS„ is used also to denote a class of 
solutions that can be transformed one to another by reflections which are the 
symmetries of equation p.2|) . We assumed also the conventions S,„ = S_,„, 
QiR^k-i) = -Q(H2fe_i), n > 0, where m is such, that < Q < 2n (0 < 
Q < 2n + 1) ioY p even (odd). To solutions that between the poles tend in 
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the limit L — oo to the vacua we will refer as pure skyrmionic solutions. 
To number eigenvalues of equation (|2.4j) we use the convention that they are 
labelled by the number of nodes of the respective eigenfunctions to which 
they corresponds, in particular the lowest eigenvalue is denoted by Aq. 



4.1 Nonsymmetric solutions possess constant number 
of modes of instability. Pure skyrmionic solutions 
are stable 

Regardless of whether a solution of ()1.2|) appears together with its compan- 
ion (as in the special case where S3 appears together with S2H1, figure El) or 
bifurcates from an already existing one (e.g Si bifurcates from HJ, this is the 
rule, that this occurs always in conjunction with vanishing of some eigenvalue 
of the Hessian. This paragraph summarizes the observations made for non- 
symmetric solutions, i.e. for solutions which have indefinite parity. There 
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Figure 3: Evolution of eigenvalues (instead of eigenvalues A here are shown 
the values A/L as functions of := 2/L^) for nonsymmetric solutions S3 

( ) and S2H1 ( ), on the left and for S1H3 ( -), S2H2 ( ) on 

the right. 



exist two cases. 

1. For a fixed triad of integers {n,k,m) a solution S„H2fe-2S„_„ appears 
together with S„_iH2fc_iS„_„. The indices are such that k > 1, n > 1 
and m is arbitrary as far as the inequality 1 < Q <2n — 2 is satisfied. 
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Then the eigenvalue X2k^2 vanishes at a characteristic for each triad 
(n, k, m) critical value L = L*^f^^ (both the solutions have the same 
spectrum only at L^km)- The i^*^^ is the same value of L at which 
this pair appears, and respectively, for S„H2fc_2S„_„ or S„_iH2fc_iS„_^ 
the eigenvalue X2k-2 positive or negative as L > L*^^. For L > L*^^^ 
the map S„H2fc_2S„_„ has 2k — 2 and S„_iH2fc_iS„_,„ has 2k — 1 modes 
of instability. The number of instabilities is the same as the number 
of instabilities of the respective harmonic maps h2fe_2 or hj^,.! to which 
these solutions tend pointwise in between the poles as L —>■ oo. In 
figure El are shown spectra of the pair S2H2 and S1H3 corresponding to 
the triad (2, 2, 2), and of the pair S3 and S2H1 corresponding to (3, 1, 3). 

2. For a fixed triad of integers (n, k,m) the solution S„H2fe_iS„_^ appears 
together with S„_iH2feS„_„. The indices are such that k > 1, n > 1 and 
m is such that the condition 1 < Q < 2n — 2 holds. Then X2k-i = 
at a characteristic for each pair its own critical value of L at which 
this pair appears. Similarly, as before, these solutions possess constant 
number of instabilities and inherit them from the respective harmonic 
maps. 

Summing up, each nonsymmetric solution possesses a constant number of 
instabilities. This number is equal to the index of a harmonic map to which 
the solution tend pointwise in the limit L ^ 00. Note also that within a 
given pair its members differ by one in the number of instabilities. There also 
exist a subclass of (P) which is composed of pairs containing pure skyrmionic 
solutions. The last are characterized by the property that they are stable, 
i.e. pure skyrmionic solutions are always stable. 

Numerics also show (and this is true for all solutions, and can be easily 
proved by examining equation ()2.4j) ) that the spectra of eigenvalues of equa- 
tion ()2.4|) with the weight function g^ip) = 47rsin^'?/', in the limit L — > 00, 
tend (if divided by L) to the spectra of eigenvalues of the harmonic maps 
between three-spheres [T]. This correspondence was the main motivation for 
the choice of this special g. Of course, there exist a more general class of 
weight functions for which such limiting behaviour would be observed (it 
should be clear that, due to the properties of our solutions, it would happen 
e.g. for the g which originates from the kinetic term ()B.2|) and from definition 
()B.4|) ). The same rescaling of energies (jl.4j) reproduces in this limit energies 
of harmonic maps [Hj . 
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4.2 Symmetric solutions S„S„, n > 1 (of parity (—1)) 

This subclass of pure skyrmionic solutions is distinguished by the fact, that 
all the solutions exist for all values of L. These solutions have parity (— 1) 
and positive definite spectra (figure E}. Due to the positive definiteness we 
conclude that the solutions S„S„ are stable for all L. Apart from the poles 




Figure 4: Evolution of shooting parameters (on the left) and of eigenvalues 
(on the right) for the stable solution SiSi. Instead of eigenvalues here are 
given the rescaled values X/L, also here is introduced the variable = 
Solutions S„S„ with n > 2 poses qualitatively similar and positive spectrum 
of eigenvalues, thus are also stable for all L 

of the base three-sphere these solutions tend pointwise to the vacua. Also for 
all n the spectra of eigenvalues tend to the universal spectrum of Hq which 
was calculated in ^3.21 

4.3 Symmetric solutions S„HiS„, n > 1 (of parity (—1)) 

These solutions exist for all L and are stable for sufficiently small L < L*^ 
where {L*^} is a monotonically increasing sequence of critical values of L. 
For a fixed n > 1 (this is also true for n = 0, i.e. & Si - the case which 
was already analysed) the transition between the stable and the unstable 
stage of the solution S„HiS„ is accompanied by the appearance of the pure 
skyrmionic solution S„+iS„, together with S„S„+i, which are not refiection 
symmetric solutions and have positive definite Hessian - thus are stable. As 
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Figure 5: Evolution of shooting parameters (on the left) and of eigenval- 
ues (on the right) for solutions SiH^Si ( -) and SjSi ( ). Instead of 

eigenvalues here are given the values \/L. Here is also introduced the vari- 
able := Solutions S„HiS„ with n > 2 possess qualitatively similar 
in behaviour spectrum, that is for some critical L the lowest eigenvalue of 
S„HiS„ becomes negative and this is associated with the bifurcation of the 
new stable solution S„+iS„. 

an example in figure El are shown characteristics of the solution SiHjSi. For 
L > L* the solution S„HiS„ inherits its single instability from the harmonic 
map hi to which, between the poles, it pointwise tends as L — > cxd. For all n 
the limiting spectrum is reproduced by eigenvalues of Hj. 

The critical behaviour of solutions S„HiS„ is quite analogous to the be- 
haviour of the identity Hi. In particular, due to the transition in instability at 
the critical L* new solutions S„+iS„ and its reflection S„S„+i, both with spon- 
taneously broken reflection symmetry, appear by bifurcation from S„HiS„. In 
more detail it is described in the next paragraph. 
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4.4 Symmetric solutions S„H2fe+iS„ of parity (-1), S„H2fcS„ 
of parity (+1), and accompanying solutions, n > 0, 
/c > 1, (thus also stable pure skyrmionic solutions 
S„S„ of parity (+1) and unstable pure harmonic so- 
lutions of parity (±1)) 

The analogous critical phenomenon, by which the solution Si was born by 
instability of Hi, we come across by observing the evolution of the reflection 
symmetric solutions S„H2fe_iS„ and S„H2fcS„. This instability is responsible 
for the appearance of new solutions by bifurcations. Unfortunately, unlike in 
the case of Si, these instabilities can not be analyzed here analytically since 
all solutions but Hi(= SqHiSo) are known only numerically. 

A solution S„H2feS„, n > 0, A; > 1 (with eigenvalues A; =: an,i) appears to- 
gether with S„+iH2fc_2S„4.i (A; =: Pn,i) at some characteristic L := L„ 2fc when 
«n,2fc-2 = = (32k-2- Both solutious are of parity +1. (If A; = 1 the last solu- 
tion is simply S„+iS„+i thus is pure skyrmionic solution). For L > L„2fc the 
solution S„+iH2fc_2S„+i has 2k — 2 modes of instability, which is the same as 
the number of instabilities of the harmonic map h2fc_2. Moreover, an,i < I3n,i 
for L > Ln,2k- For Ln,2k < L < L*^ 2k the eigenvalue an,2k-i is positive, 
tends to as L — > L* 2k and is negative for L > L*^ 2k- Thus, for L > L*^2k 
S„H2fcS„ has 2k modes of instability - the same as the number of modes of 
the harmonic map hafe. At the critical L* 2k the solution S„+iH2fc_iS„ and, due 
to reflection symmetry, the solution S„H2fc_iS„+i appear by bifurcation from 
the solution S„H2fcS„. As an example in figure IHl are presented parameters 
of the solutions Ha, SjSi and SjHi. The solutions S„H2fcS„, > 0, > 1 
and S„+iH2ft_2S„+i differ by 2 in the number of negative eigenvalues if L is 
sufficiently large. This is the first solution which possesses 2 modes of insta- 
bilities more than the second solution, and which 'bears' the new solution. 
The new solution has 2k — 1 modes of instability (the same as the map h2fc_i) 
and is nonsymmetric, unlike the one from which it bifurcates. 

The solution S„H2fc+iS„, with eigenvalues A; =: a'^i, appears together 
with S„+iH2fe_iS„+i, whose eigenvalues are A; := at some characteristic 
L = Ln,2k+i at which a'n2k-i = = P'n2k-i (figure [7j). These solutions have 
also definite parity, i.e. (—1). For L > L„ 2A:+i the solution S„+iH2fc_iS„+i has 
2k — 1 modes of instability, which is the same as the number of instabilities of 
the harmonic map h2fe_i. Moreover ; < ^ for L > Ln,2k+i- For Ln,2k+i < 
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Figure 6: The evolution of shooting parameters and eigenvahics for sohitions 

Ha ( -), SiSi ( ), and SiHi (with its partner HiSi ) ( -). This 

is a particular case of more general situation when S„H2feS„ (n > 0, /c > 1), 
together with S^+iHjfc.aS^+i appear and, due to the transition in the number 
of instabilities of S^Hj^-S^, the solution S„+iH2fc_iS„ appears by the process of 
bifurcation. All the solutions have vanishing topological charge. 

L < L* 2k+i the eigenvalue 2k positive, tends to as L ^ L* 2k+i ^-nd 
is negative for L > L^2k+i- Thus, for L > L'^2k+i solution S„H2fc+iS„ 
has 2k + 1 modes of instability - the same as the number of instabihties 
of the harmonic map hj^+i. This is the critical value -^^ 2fe+i which the 
solution S„+iH2fcS„ (and S„H2fcS„_,_i) appear by bifurcation from S„H2fc+iS„, 
i.e. from this one of the original pair, which has 2 more, than the second, 
modes of instability. The new solutions have constant number of 2k modes 
of instability (the same as ha^,) and have spontaneously broken parity. 

5 The mechanism of the spontaneous break- 
ing of the reflection symmetry 

The instabihties described above are closely related to the breaking of re- 
flection symmetry. As the rule we state that solutions which bifurcate from 
symmetric solutions are not symmetric and that only the solutions which 
are symmetric turn out 'to be able to give birth' to new solutions. The 
special class of hedgehog solutions is described by the reflection symmet- 



26 




Figure 7: Evolution of eigenvalues for solutions H3 ( -), SiHiSi ( ), 

and S1H2 (H2S1) ( -) on the left, and to compare with, the same character- 
istics for S1H3S1 ( ), S2H1S2 ( ), and S2H2S1. Instead of eigenvalues 

A, are shown A/L, here is used also the variable k"^ = 2/L?'. Both figures are 
misleadingly similar, note the difference between the critical values of and 
tiny differences between the critical values of eigenvalues in both figures. Be- 
cause of this one would have the impression that the presence of skyrmions 
has not very crucial influence on the stability, but it suffices to compare this 
with the spectrum of S1H3 (figure Ej) - which is not a symmetric solution - to 
see this is not the case. These solutions are the particular case of the more 
general situation when S„H2fc+iS„ {n > k > 1) together with S„+iH2fc_iS„+i 
appear and, due to the transition in the number of instabilities of S„H2fc+iS„ 
the solutions S„+iH2fcS„ (and S^HjfeS^+i) appear by the process of bifurcation. 
These solutions are topologically nontrivial Q = 2n + 1. 

ric equation ()1.2p . but some solutions possess no definite parity. This phe- 
nomenon is known as the spontaneous breaking of reflection symmetry or of 
parity. In the case of the Skyrme model on the following reflection sym- 
metries of equation ()1.2|) are being spontaneously broken: F{x) F{—x) 
and F{x) Qtt — F{—x) (we assume lim^^_oo F{x) = 0). We used here the 
variable x = Intan ('?/'/2). The second symmetry is a composition of several 
primitive symmetries. Because of this we ascribe a definite parity +1 or —1 
to a solution which is invariant with respect to the first or with respect to 
the second transformation. By examining the process of the appearance of Si 
by bifurcation from Hi, which was performed analytically, the phenomenon 
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of spontaneous breaking of the reflection symmetry can be understood. By 
comparison with numerical data of other solutions we can generalize our ob- 
servations. As it will be shown, the explanation of the symmetry breaking 
can be successfully carried out with the use of the modes whose eigenval- 
ues vanish at some critical radii. This construction follows straightforwardly 
from the analytical approach of ^3.31 Here it should be noted that the first 
time the wording 'spontaneous breaking of full rotational symmetry' in the 
case of Skyrme model on was used, was in the paper P . 

By a critical mode we will call the mode to which it correspond the eigen- 
value which vanishes at some critical L = L*. We observed that for L < L* 
this eigenvalue is positive and negative for L > L* in all cases. It is worth 
to note the fact that a critical mode has always definite parity at least for 
L = L*. This is because the solutions that 'give birth' to other solutions 
by bifurcations have definite parity, and we observed that only such solu- 
tions possess critical modes. The symmetry of the Sturm-Liouville equation 
()2.4j) . which determines these modes for such solutions, can be broken only 
by the weight function which can be arbitrary. In fact, the vanishing of an 
eigenvalue means that the corresponding eigenfunction is not affected by the 
weight function, and by this is intimately connected with the Hessian alone 
and consequently with the equation ()1.2|) . Thus, by analogy with ^3.31 and 
by the observation of similar behaviour of shooting parameters in the vicin- 
ity of critical points i^* p, at which new solutions appear by bifurcations, 
we can write down the formulas that describe, in the first approximation, 
the shooting parameters of solutions S„+iH2fe_iS„ and S„H2fc_iS„+i (which 
bifurcate from S^Ha^S^) or S„+iH2feS„ and S„H2fcS„+i (which bifurcate from 
S„H2fc+iS„). For sufficiently small L > L* ^ these formulas read respectively 

S„H2fcS„ ± 7„,2fc^n,2fc-i or S„H2fe+iS„ ± 'jn,2k+iCn,2k whcrc ^n,p IS the critical 
mode and the coefficient 7„_p shows the characteristic for critical phenomena 
rational power dependence on (L — L^p) (the exponent was calculated in 
S (|3.3p to be exactly 1/2 and with very good accuracy this value was also 
derived by fitting curves to numerical data for several solutions). Note also, 
that the mode ^n,2fc-i has an odd number of nodal points (hence the parity 
— 1) while the solution S„H2ftS„ has parity +1. Analogously ^n,2k has parity 
+1 which is opposite to the parity of S„H2fe+iS„. Thus, the solutions which 
appear by bifurcation process have indefinite parity, i.e. have spontaneously 
broken refiection symmetry. In the first approximation, it is due to the ex- 
citation of the critical mode which has the opposite parity with respect to 
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the parity of the excited reflection symmetric solution. This occurs while 
the qualitative transition takes place, when the reflection symmetric solution 
undergoes increase by 1 in its number of instabilities. In other words, one can 
say that at the critical point, in a sense, it costs no energy of the excitation 
of the critical mode. This is always inevitably connected with the process of 
spontaneous breaking of reflection symmetry. 

6 Conclusions 

In this paper the linear stability analysis of the whole spectrum of the 
Skyrme's hedgehogs on the three-sphere was carried out. We considered only 
spherically symmetric perturbations. The most remarkable it proved the 
influence of instability modes on the existence of a class of solutions with 
spontaneously broken parity and which appeared by a critical phenomenon. 
Also some general remarks on the construction of the Hilbert space of admis- 
sible perturbations were presented together with the resulting problem of a 
choice of a norm and its influence on the resulting spectrum of perturbations. 

The main lesson of this paper is that for L sufficiently large (L is inter- 
preted as the radius of the base three-sphere) the number of instabilities of 
solutions S„HfeS,„ is the same as the index of the harmonic map h^. to which 
(apart from the poles) the solutions tend pointwise as the radius grows. The 
appearance of new solutions is always accompanied by the vanishing of some 
eigenvalue of the spectrum of the Hessian. In particular, pure skyrmionic 
solutions S„S„ localized at the poles are always linearly stable. The num- 
ber of instabilities of nonsymmetric solutions is independent on L. There 
also exist a class of reflection symmetric solutions whose number of negative 
modes increase by one as L grows. Due to the qualitative transition new 
solutions with broken reflection symmetry appear by bifurcations and this 
process possesses many characteristics of a critical phenomenon. In partic- 
ular, the identity solution H^, which is stable for small L, becomes unstable 
at L = \/2. Then the 1-skyrmion Si (together with its reflections) is born by 
biffurcation from the Hi and exists as a stable solution for L > \/2. In this 
unique case the critical phenomenon may be astonishingly simply described 
analytically and therefore completely understood. This gives rise to a general 
picture of the process of the spontaneous breaking of the reflection symme- 
try which may be intuitively presented as a costing no energy inflation of the 
critical mode on the unstable symmetric solution. The unstable solution and 
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the critical mode have mutually opposite parity, thus the newly born solution 
has spontaneously broken parity. Finnaly, we also found two unique series 
expansions for the profile and for the energy of the 1-skyrmion (although the 
respective expansion coefficients are not known in general forms) which read 
respectively 

, , 2 3 , 29369 

F{f)=f + xsmf + x—sm2f-x I g^gggg 

4 f 4217 , 109 , \ 

—X sm 2w sm 4w + 

V887040 ^ 20160 

and 

E\F] 3^2 / 11 2 209 4 5209 « 2126641283 o ^ 

— - — - = — — I 1 H X X H X X + 0(x 

127r2 4 V 180 10800 864000 1149603840000 ^ 

where 

These series expansions are valid in some right neighbourhood of the critical 
radius L = \pl and numerical assessment of their behaviour makes them 
plausible to possesss nonzero radii of convergence. 

This critical behaviour which we observed in our model exist neither in 
the Skyrme model on the fiat space nor in the model of harmonic maps be- 
tween three-spheres when free parameters of these models - the characteristic 
size of a soliton in the first case, or the radius of the base three-sphere in 
the second case - are being changed. It shows that by coupling together two 
field theories we may produce another theory in which new phenomena may 
appear that are absent in the decoupled case. Dimensional parameters intro- 
duced by different models determine then some dimensionless numbers whose 
values are usually crucial to the existence of different types of solutions. This 
fact is known when nonlinear fields are coupled to gravity 0^ . The main mo- 
tivation of our work was to understand such model-independent phenomena 
in possibly the simplest case which is a nonlinear scalar field theory on a 
fixed space-time background. 
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A Formal series expansion for the profile and 
energy of the skyrmion Si on 

In order to find a formal series expansion for the profile F{tp) of Si in the 
vicinity of L = \/2 we assume that this solution can be written using the 
following ansatz: 



F{'^,x)=^ (A 

2 / , sin3^/'\ A f-, sinS?/' sin5^/'\ 

ai + X 03 1 + 03,3^ — - + X 05 1 + 05,3^ — - + 05,5^ — - + . . . 

sm ip J \ sm ip simp 

■X sin 

, 2 / sin4'?/'\ 4 / sin Alp sm6ip\ 

+ <^ 02 + X 04 1 + 04,4 . + X ag 1 + a6,4 . ^ , + ae,6 . ^ , + 

\ sm zip J \ sm zip sm zip ^ 

„2 



X sin 2-0 



where 



X 




L > V2, and ai := 1. 



This ansatz follows from the idea that the solution F{iIj, x) should be linearly 
decomposed in the base of eigenfunctions of the Hessian (j2.ip evaluated at Hi. 
This is implied by the appearance of Si at the critical radius L = \pi by the 
excitation of the lowest energy mode (i.e. sm.i\}) on Hi, which was perfectly 
confirmed when we reproduced analytically the numerical observation that 



30 V2 / 

ia-i| = y 

which implies that 
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L-^Pl\■ 



X, 



dF{ijj, x) 



dx 



= sinip, 

x=0 



and the vague idea that this behaviour should be somehow analytically con- 
tinued in the variable x. From the symmetry of equation p.2|) we impose 
on the series also the requirement that, since L oc x^, the second solution 
7T — F{ti — Ip, x) should be simply F{il), — x) (both skyrmions bifurcate from 
Hi at X = and can be transformed one to another by the reflection with 
respect to the point (7r/2, 7r/2)). 
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The coefficients and ; (at least the one we had found) are rational 
numbers and relations between them become more and more complicated 
as k increases. Due to nonlinearity of equation ()1.2|) it seems there is no 
general recurrence which would determine these coefficients. Anyway, the 
above ansatz enables us to find them all successively step by step. In table 
n we give all the coefficients which are required to determine the profile 
of Si with the accuracy up to 0{x^'^). The question arises if the formal 
series possesses nonzero radius of convergence. We gave in the main text 
some arguments that this series is really a solution. In order to enhance 
our argumentation we may try to estimate this radius. We can treat the 
maximal absolute values of the functional coefficients standing in braces at 
in equation ()A.1|) to built a majorizing power series in the variable x. 



X 



2k 



Let denote the respective values by {ci, C3, . . . } and {c2, C4, . . . } respectively 
for the first and the second brace. Thus the partial sums of the above series 
are majorized by ^c„|a:;|". In figure |H1 are shown 20 initial elements of the 
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Figure 8: 

left: Initial elements of the sequence {1/ ^yc^}, where X^n*^"-^"^ ^ majoriz- 
ing series for the profile of Si. This strongly suggest that the Taylor- Fourier 
series expansion ()A.1|) has a nonzero radius of convergence and even that, it 
is absolutely convergent. 

right: Initial elements of the sequence {1/ V^2n} where are defined in 
flA.2|) . This suggest that the Taylor series expansion of energy of Si may have 
a nonzero radius of convergence. 



sequence (1/ -^/c^). Their behaviour strongly suggest the hypothesis that the 
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Taylor-Fourier series expansion ()A.1|) has a nonzero radius of convergence 
and even that, it is absolutely convergent. 

The Taylor series expansion of the energy of Si around x = must con- 
tain only even powers of x since F{ip, x) and F{ijj, —x) must have equal 
energies. Therefore the energy of Si is a smooth function of L — a/2, unlike 
the F'(0,x(L)) — 1, which has a branch singularity at L = v^. This series 



and the respective coefficients are given in table |21 For x = we get the 
energy of Hi at L = \/2. Here we also construct the sequence {1/ V'^2n} to 
make it plausible that the series ()A.2|) and, indirectly, the series ()A.1|) . have 
nonzero radius of convergence. In figure |H1 are shown initial elements of the 
sequence (1/ V^2n)- Their behaviour suggest the hypothesis that the Taylor 
series expansion ()A.2|) may have a nonzero radius of convergence. 

B On the weight function originating from a 
kinetic term 

In this appendix we derive the formula for the weight function, which natu- 
rally follows from the calculus of small oscillations about equilibrium. We ob- 
tain this by the comparison of two conceptually different approaches to linear 
stability analysis. Suppose, for example, that we have a differential equation 
for a function of one space variable and that the equation is a condition for 
extrema of some functional. If the function is a solution and, moreover, a 
local minimum of the functional then it is also a stable equilibrium of some 
1-1-1 dimensional field theory. Stability is required for physical reasons. In 
order to check that the function is a minimum it suffices to prove positive 
definiteness of the second variation of the functional. This in turn requires 
the existence of a scalar product (which should be defined by additional ar- 
gumentation) and defines a Sturm-Liouville problem. On the other hand, in 
linear approximation, the calculus of small oscillations about the equilibrium 
produces another Sturm-Liouville problem with a definite weight function 
determined by the kinetic term of the field theory. It turns out that both the 
problems comprise identical self-adjoint differential operators (the function 
spaces of admissible perturbations are the same by construction). Thus, to 
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ai = 1 



«2 - +m 



_ 29369 
3 316800 

7260 



^3,3 



29369 



_ 4217 
^ 887040 

4796 
4217 



^4,4 



as = + 



12335768813 
1405071360000 

_ _ 4404867060 
"5,3 — 12335768813 

3226815900 



_ _ I 5010043861 
"6 — "T 7376624640000 



06,4 



33475541308 
5010043861 



05,5 



+ 



12335768813 



_ _ I 13362918261 
"6,6 — -r 10020087722 



_ I 7986485836639 
"7 "'"25751126016000000 

158986289932080 
7986485836639 

57715917199500 



+ 

"7,3 = + 
07,5 
"7,7 



+ 



7986485836639 

6755487649680 
7986485836639 



_ 473077794653026103 
"8 98150416809984000000 

_ 361662885996135944 
473077794653026103 
_ I 51442637682491484 
"8,6 "I" 473077794653026103 

_ 5773312055063520 
"8.8 473077794653026103 



"9 = 

"9,3 

"9,5 

"9,7 

"9,9 



+ 



127577656118454895711 
167510044689039360000000 

790229840916734772664 



+ 



127577656118454895711 

309555559736002696600 
382732968355364687133 

20023176097891467424 
382732968355364687133 

2544822818150368800 
127577656118454895711 



_ I 4211381544228526601549071 
"10 "•" 1061084526049375027200000000 

_ _ 822641146563016022074558 
"10,4 — 4211381544228526601549071 

_ 1447811040980803898316621 
"10,6 33691052353828212812392568 

_ 13895432790539792801760 
"10,8 4211381544228526601549071 

_ I 55236971873923883695575 
"10,10 — "T 33691052353828212812392568 



_ 107845172239555936226487785603477 
"12 73415504602973334681944064000000000 

_ I 60201002753190711583565881906936 
"12,4 -r 107845172239555936226487785603477 

„ _ 8388700062249632787593687120376 
"12,6 — 107845172239555936226487785603477 

„ _ 5165826035166448141662612158400 
"12,8 107845172239555936226487785603477 

_ I 975465355644891981999033888060 
"12,10 -r 107845172239555936226487785603477 

_ 68171833635114106739516947272 
"12,12 107845172239555936226487785603477 



"11 = 
"11,3 
"11,5 ■ 
"11,7 
"11,9 ■ 
"11,11 



+ 



982955548821103176200169976183 
2390411220284032061276160000000000 

I 3024748425995528913667880277780 



982955548821103176200169976183 
1248755928443601042244324285400 
' 982955548821103176200169976183 

58074313732006509362657835680 
982955548821103176200169976183 

56310543306011901706414401600 
982955548821103176200169976183 



+ 



68232461 16254206530384345600 
982955548821103176200169976183 



_ 1766865169954398109033856353883942537617 
"14 6132984423523186432660243218432000000000000 

_ 4349592583735535844513687814933733037976 
"14,4 1766865169954398109033856353883942537617 

_ 1508853846688227478046216373401843325339 
"14,6 1766865169954398109033856353883942537617 

_ I 585277715345242542258046459394514314520 
"14,8 — -I- 1766865169954398109033856353883942537617 

_ 78074347310444368596245949472056202125 
"14,10 3533730339908796218067712707767885075234 

_ I 1726765295922525300215513807655077490 
"14,12 -I- 1766865169954398109033856353883942537617 

„ _ 38714472844152751765544208844685625 

"14,14 199083399431481477074237335648894933816 



"13 = 
"13,3 = 
"13,5 = 
"13,7 = 
"13,9 = 
"13,11 
"13,13 



43153404274707851225466988047569832167 
413476121923945820928708968448000000000000 

313572180766907863870688164387239391620 



+ 



43153404274707851225466988047569832167 

137371166634494138333962660336190541100 
43153404274707851225466988047569832167 

65940526790056736560082171865362030240 
43153404274707851225466988047569832167 

16887408629502465138357459627899436000 
43153404274707851225466988047569832167 

1437225091006557895313712800717160000 



43153404274707851225466988047569832167 

92252259291662946944254093579272000 
'43153404274707851225466988047569^|2167 



Table 1: Initial coefficients of the Taylor- Fourier series expansion ()A.1|) of 
the solution Si in the vicinity of the critical radius L = a/2- 
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"6 "T 864000 

^ _ 2126641283 
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1149603840000 

_ , 8434513426021 
10 — "T 16995743170560000 



_ 464402526814953371 

12 — 6730314295541760000000 

_ 23242947675020983465203089 

14 — 626856089235323092992000000000 



J _ I 14367315331877365555523427348253 
"16 — +- 



il6 -r 379122562769523406641561600000000000 



31939951417910304147919969132085999 



2017781268474468255763984809984000000000 
' _ 74250563114418365106618380718783299763796141 
20 109069148686118907097066434918875136000000000000000 



Table 2: Initial coefficients of the Taylor series expansion of energy of Si at 
the critical radius L = ^J2. In order to calculate d^^ we need the proffie of 
Si to be determined with the accuracy up to x^". 



specify the weight function in the ffist case, it is natural the requirement that 
both methods should give the same results. Of course, to check positive defi- 
niteness of the second variation of the functional we can assume an arbitrary 
weight function. If the function is not a minimum we are not guaranteed 
that the different field theories have qualitatively the same dynamics in the 
vicinity of the common equilibrium. It would be so, under the condition 
that the number of instable modes were independent of the weight functions, 
which is not obvious. In what follows we determine the weight function. 

We can supplement functional ()1.4|1 with an arbitrary kinetic term not 
depending explicitly on time (to have the total energy £^[F(r)] conserved) 

e\F{r)\= j wii;) (^K{^,Fii;,T),Fii;,T)) + V{^,Fii;,T),F'{^,T)))dij, 



and observe the resulting time-evolution of some initial data in time r. Con- 
ventionally we measure duration using natural units of length e~^f~^ then 
r is dimensionless. The quantities K and V are interpreted as densities 
of kinetic and potential energy respectively. We assume that K = KF^ 
where K is positive function of , F, F' , F . A perturbation can be imagined 
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as a time dependent excess from the equilibrium. An equilibrium solution 
F is said to be dynamically stable if it remains within some finite bounds 
||5F(r)|| < £:(A) for every r if only the initial perturbation is sufficiently 
small such that 5E < A. Thus, to check stability an initial solution is re- 
quired to be not arbitrary but such, as its total energy was slightly above 
the energy of the equilibrium solution and in its pointwise vicinity. It is 
clear that different kinetic terms give rise to different time evolution and the 
equilibria are common and the same as solutions of equation ()1.2|) . since the 
condition that 5U[F]1 5F{ip) = for static solutions is the straightforward 
consequence of equation of motion 

dK_d_dK_6U[F]^Q 

dF dtdF SF ^ ■ ^ 

and of the requirement that the solutions were in equilibrium. If we had 
decided to assume the kinetic energy density as in the original Skyrme model 
on S^, i.e. 

K{^,F,F)=(l + ^'^)l^F\ (B.2) 
\ L sm ip J 

then the evolution of the system would have been governed by the equation 

(^L sin^ ^ + I sin^ F^ {f" - L^F^ + sin 2F {f'"^ - L'^F^^ + L sin 2ijF' 

1 sm^F\ 
L + -—^ sin2F = 0, 
L sm ip J 

while assuming K{iIj, F, F) = L^F"^ we would have been led to the equation 



L sin^ + - sin^ F ] F" - sin^ tpL^F + (L + F' sin 2z^F' 

L J \ L sm2yj J 

( Isin^FA 

\ L sm '^} J 

As it will turn out the last choice for K gives, up to a constant factor L^, the 
natural on weight function 47rsin^?/'. We can always look for solutions 
in the form of F{tp,T) = F{iIj) + eri{il),T) where e measures the scale of 
perturbation, i.e. its energy, which is proportional to e^. Next, for small e by 
analogy with the calculus of small oscillations, solutions (if exist) for which 
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rj is bounded for all r, may be arbitrarily well approximated by the following 
linear equation, which is derived by equating with zero the term of the order 
0{e) in the Taylor series expansion of equation ()B.1|) in the variable e, i.e. 

d^U dri\' f d^U f d^U \'\ d^Kd^T] , 

""dF^d^ ) + r ^ ~ VdFT J J ^ = -""oFw ^ ^ ^ 

Other terms (not shown) which enter above equation vanish since we per- 
turb a static solution for which F = 0. If the equilibrium was not stable the 
approximation would still have been useful to predict the existence of insta- 
bilities. It follows, that for the class of kinetic terms assumed it is always 
possible to find r] by the method of separation of variables rj{ip,T) = ^('(/')a(r), 
where a is of the form a(r) = e*'^'^ and cu^ < for unstable modes. Thus 
the dynamics impose the natural choice that weight function g{ip) should be 
defined according to the formula 



giijj) = wiip)-- 



2dFdF 



(B.4) 



equil 



since equation ()B.4|) reduces to equation ()2.3p after the substitution rj{ip, r) = 

sin {uT + 0) with A = o;^ or, if the mode is unstable, rj{ip, r) = ^{ip) ^xp (icur) 
with A = —uj^. Put differently, multiplying ()B.3j) by 77/2 and integrating, 
yields 



2 ^ J^^J J 2dFdF^ 



2 

since ci + Aa = which reproduces ()2.3j) for the special choice of g. In 
this sense, the spectrum of fj2.4j) with g = w may be interpreted as 0;^ for 
a spherically symmetric scalar field on in an equilibrium if = F"^. 
Thus, the Hilbert space spanned on vanishing on boundaries and mutually 
orthogonal (with respect to the special g) eigenf unctions of the self-adjoint 
equation ()B.3j) is constructed and each perturbation r/ is a superposition 
r]{iJj,T) = ^^^kWiakakir) + hkl3k{T)) where denotes sinwfcT or sinhwfcT 
and [3k denotes cosa;fcr or coshcUfcT, respectively, if A^ is positive or nega- 
tive. Only a finite number of modes can be unstable. The total energy of 
the perturbation is finite and quadratic in amplitudes (proportional to e^) 
and, of course, only if all A^ are positive the function F{;ip) + eri{ip,T) well 
approximates the particular solution of equation ()B.1|) for all r. 
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It is thus seen that the hnear dynamical stabihty is closely related to 
the (abstract) notion of linear energetical stability defined in 12.2.11 and, if 
the weight function g is appropriately chosen, all eigenvalues of the Hessian 
may be directly interpreted as squared eigenfrequencies ujI of a 
system whose dynamics is related to some kinetic term. 
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